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TECHNICAL NOTE 4278 

APPLICATION OF STATISTICAL THEORY TO BEAM-RIDER 
GUIDANCE IN THE PRESENCE OF NOISE. 

II - MODIFIED WIENER FILTER THEORY 1 
By Elwood C. Stewart 

SUMMARY 


A study has heen made of the application of Newton's modification 
of the Wiener filter theory to the optimization of a beam-rider guidance 
system operating in the presence of glint noise. Target and missile 
motions are assumed to be coplanar. 

The theory has been applied using a typical variable-incidence 
missile by placing a realistic restriction on the mean-square surface 
deflection so that the system operation is confined to the linear range. 
The transfer functions of the optimum guidance system are derived and 
an example missile-control system is synthesized. It is shewn that the 
minimum attainable error corresponding to a realistic control restriction 
is close to that for the Wiener theory. Performance of the system versus 
noise magnitude is given. 

It is shown that the most critical saturable quantity is control 
deflection since constraining control deflections to realistic values can 
prevent saturation of other important quantities. Servo energy require- 
ments are also greatly reduced in comparison with systems in which 
saturation is allowed to occur. 


INTRODUCTION 


Noise effects in missile-guidance systems impose one of the most 
serious limitations on the effectiveness of a missile. Noise signals in 
general have the effect of adding to the system false information which 
cannot be distinguished from true information. Thus the missile responds 
to the noise signals as well as to the true signal; the miss distance is 
thereby increased. Since many sources of noise (such as angular scintil- 
lation) are inherent in the physical mode of operation and cannot be 
removed, it is important that missile systems be designed to minimize the 
miss even though the noise is present. Statistical filter theory is useful 
in performing this minimization. 

One of the most valuable theories has been developed by Wiener 
(refs, 1 and 2). A previous report has considered the application of this 
Supersedes recently declassified NACA Research Memorandum A55Ella 
by Elwood C. Stewart, 1955 - _ 
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theory to the minimization of the effects of radar glint noise in a beam- 
rider guidance system (ref. 3)* It was found there that the optimum per- 
formance specified by the theory could not be obtained because of certain 
practical restrictions, as for example, limiting of control deflection. 
Further study showed, however, that even in the presence of limiting, the 
optimum performance could be approached by the addition of network filter- 
ing to the guidance system. On the other hand, the optimum form for the 
added filtering could not be determined from this study since the system 
design was based on a theory which did not take into account the finite 
range of operation of the saturable quantities. In the present report an 
attempt is made to overcome this deficiency by determining the optimum 
transfer function under this last restriction. The theory for such an 
approach was made available in a paper by G. C. Newton (ref. 4) on a 
modification of the Wiener theory. 

In the application of this theory to the missile guidance problem it 
is necessary to make certain as sumptions . The main assumption made in 
reference 3> that is, the assumption that the target and missile move in 
the same plane, also will be made here. Other assumptions, such as those 
relating to the class of target maneuvers and noise, are discussed in the 


SYMBOLS 

transfer function of optimum compensating network 
transfer function of the fixed network 

noise magnitude or zero frequency spectral density, f t 2 /radian/sec 
time constant of the noise spectrum shaping filter, sec 
optimum closed-loop transfer function 
acceleration of target maneuver, ft/sec 2 

twice the average switching rate of target acceleration, l/sec 
target displacement from a space reference, ft 

apparent target displacement from true target center due to noise, 
ft 

missile displacement from a space reference, ft 

cross-spectral density of the input "signal with the desired output 
spectral density of the input signal 

spectral density of noise displacement y^, ft 2 /radian/sec 
spectral density of target displacement y T , f t 2 /radian/sec 
error between target and missile position, y T - y^, ft 
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€rp component of error e due to target motion, ft 

e-jj component of error e due to noise, ft 

p Lagrangian multiplier 

5 control- surface deflection, radians 

ljr angle of yaw, radians 

l_i 0 optimum open-loop transfer function 

[i open-loop transfer function of system approximation to (i Q 

oj angular frequency, radians/sec 

( ) complex conjugate of ( ) 

GENERAL CONS 3DERATIONS 


Glint noise is a term that is used to describe a shift in the appar- 
ent target location as determined by a radar. It is due basically to the 
variable reflection characteristics of aircraft targets and arises from 
the relative movement of the various reflecting surfaces. Since the radar 
utilizes the reflected signal to determine target location, variations in 
the reflected signal are interpreted by the radar as shifts in the target 
location. Of the many sources of noise which may exist in a missile- 
control system, glint noise is one of the most serious since it has a 
physical origin which cannot be eliminated, imposing a fundamental limi- 
tation on missile effectiveness. The situation in regard to glint noise 
is illustrated in figure l(a) where the true target position is indicated 
as yij and the glint noise is represented by the displacement y^. 

The present report is restricted, for the sake of simplicity, to a 
two-dimensional study in which the target and missile move in a horizon- 
tal plane as shown in figure 1(a). 2 The attack situation, that is the 
angle between missile and target velocity vectors, may be arbitrary. How- 
ever, in order to illustrate the theory, the tail chase illustrated in 
figure l(a) will be considered. It should be noted that for this case 
the reference line is fixed in space. The guidance system which is con- 
sidered to be of the beam-rider type is illustrated in figures l(a) and 
l(b). Its function is to make the missile position coincide as closely 
aa possible with the true target position y^. The difference y T - y M 
is indicated on this figure by the error €, which obviously should be 
minimized. The criterion of merit which will be used here is the conven- 
tional mean-square time average of the error. This criterion is particu- 
larly appropriate in the case of the beam-rider system because the target- 
to-missile range is not normally transmitted to the missile; since the 
missile never knows when the target will be reached, it is reasonable to 
mininrlze the error for all values of time. 

2 The complete three-dimensional problem would require a more complex 
analysis than used herein. Possibly either the present theory or Wiener’s 
theory for multiple time series (ref, l) could be applied to this case. 
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Emphasis should he made that although the example case is taken to 
be a tail chase, it contains the essential ingredients of the problem. 
For other attack situations a coordinate system can be defined such that 
displacements and accelerations are measured from a new reference line 
which is not fixed in space. This line is one which translates with con- 
stant velocity and without rotation. There are several suitable condi- 
tions under which such a reference line can be so defined. In the first, 
wherein the launching aircraft flies a collision course with the target, 
the line of sight between launcher and target (which coincides with the 
radar beam, ideally) would be a suitable reference line. Another case 
would be one in which the tracking radar is sufficiently far from the 
target that rotations of the line of sight are small. In either case, 
all displacement quantities as well as accelerations are then measured 
perpendicular to the reference line just as shown in figure 1(a). 


The design of a system normally depends on the inputs to be encoun- 
tered, in this case, the target motion and the noise. Because of their 
random nature, neither of these quantities can be conveniently described 
explicitly as functions of time, and statistical descriptions are more 
suitable. In the following paragraphs the target motion and noise are 
briefly described. - — . . .. 


The glint noise can be defined by specifying both the power spectrum, 
or power spectral density as a function of frequency, and the amplitude 
distribution. Many measurements have been made on the characteristics of 
glint noise. (For a brief bibliography see ref. 3 .) Although these meas- 
urements are somewhat complicated and uncertain, it is generally found 
that the amplitude distribution is approximately Gaussian and that the 
spectral density can be adequately represented by 


% = 


N 

% 2 w 2 + 1 


( 1 ) 


An examination of glint noise spectra indicates that the break point 
(l/2rtTN) is generally on the order of several cycles per second and as 
in reference 3 will be taken to be 6 cps, corresponding to Tn = 0.0265 
second. The magnitude of the spectrum, N, of equation (l) depends on 
factors such as target size and target aspect so that the guidance system 
is generally forced to operate over a wide range of magnitudes. This 
range may extend from 7 f t 2 /radian/sec for small targets up to around 
30 ft 2 / radian/sec for large bombers. For this range of magnitudes it 
might be thought that it would be necessary to optimize the system for 
each noise magnitude. As shown in reference 3, however, this was not the 
case for the method used therein since near optimum results were obtained 
by optimizing only for a raid- range value called the design value. It is 
not unreasonable to believe that the same will hold for the method of the 
present report. As in reference 3 > a. design value of 15 ft 2 /radian/sec 
will be used here. 


The type of target maneuver upon which to base the system design can 
never be determined with certainty, since the target quite obviously may 
maneuver in many different ways. A reasonable situation might be one in 
which the target pilot is only aware of the attack and therefore maneuvers 
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in some random manner to avoid being hit. One of the most useful concepts 
in system design and the one which is used in this report is to picture the 
target evasive maneuver as a stationary random process in which the target 
turns at its maximum possible rate alternately in opposite directions with- 
out regard to what the attacking missile is doing. The length of time 
between switches is assumed to be randomly distributed accor din g to a 
Poisson distribution [(l/T)exp(-T/T)] , where T represents the length of 
the interval and T the average interval length. It can be shewn that 
the spectral density of the target accleration is described by 

. = ka T 2 

^T rt(to 2 +k 2 ) 


Here the quantity arp represents the magnitude of the target accelera- 
tion normal to the beam, and k is twice the average switching rate of 
k = 2 /T. The spectral density of the target displacement is then given by 


®T = 



( 2 ) 


For the example tail- chase situation, the target is assumed to maneuver 
with ±lg acceleration with an average period of five seconds, which gives 
Qjjt = 32.2 ft/sec and k = 0.4 switch/sec. 

There are several important comments to be made concerning this type 
of input: (l) First of all, a statistical description of the target maneu- 

ver process is a desirable one, since target motions cannot be described 
as unique functions of time. (2) It is clear that the maneuver assumed is 
a severe one and puts the system to a good test; it is often found that 
systems designed according to theories based on either no maneuver or very 
weak maneuvers have unacceptably poor performance in the presence of a 
more severe maneuver. (3) Another consideration not generally realized is 
that the stationary process described above is also applicable to certain 
important nonstationary processes. In any real problem it is apparent 
that the inputs are distinctly nonstationary. For instance, they are non- 
stationary because the target motion and noise do not exist for an infi- 
nitely long time into the past. However, the nonstationary character of 
the input is due to the strict mathematical definition. It is clear that 
in the practical case it makes little difference to the missile, so far as 
miss distance is concerned, whether a process persists over an infinite or 
a finite period so long as the process begins at a time before the end of 
the attack by an amount equal to or greater than the system response time. 
(Of course, the process may terminate any time after the attack is over 
without affecting the results.) In other words, an infinite period is, 
for practical purposes, simply one which is longer than the system response 
time. Thus, when the system response times are short, results obtained by 
means of the stationary input apply directly to an important class of non- 
stationary problems. The results presented herein are in this category. 

(4) It should be pointed out that an optimum design based on the above tar- 
get maneuver will operate efficiently against this class of maneuver as a 
whole, although not necessarily efficiently against any one particular tar- 
get maneuver (e.g., a single target turn). However, as indicated in refer- 
ence 3, this design will be essentially optimum for the single- turn maneuver 
as well. 



6 


NACA TN 4278 


ANALYSIS , RESULTS, AND DISCUSSION 
The Modified Wiener Theory 


The primary objective of the modified Wiener theory is to arrive at 
an optimum transfer function in which the effects of certain saturating 
elements are considered. Ideally it would be desirable to determine this 
optimum transfer function without any assumptions about whether or not the 
system would operate in these saturation regions. However, no such theory 
is available at present. The approach considered herein is one in which 
the transfer functions of certain available elements are chosen so that 
saturation does not occur. Thus the system is forced to remain linear 
and linear methods are immediately applicable. Although It iB not appar- 
ent at first glance, there are certain indications that this approach may 
lead to desirable results. First, it appears intuitively that saturation 
is undesirable because of the loss of both intelligence and ability to 
control the missile when in the saturated regions. Second, it iB indi- 
cated In reference 3 that the undesirable effects of saturation can be 
reduced by the addition of filtering without undue increase in the error; 
this method in effect tends to keep the operation of the system in the 
linear range. These are only indications, however, and the real justi- 
fication in the completely linear approach depends on the resultant error 
performance which can only be established by investigation. 

The theory to be used in this approach follows Newton’s modification 
of the Wiener filter theory. To describe the essentials of the theory it 
will be convenient to refer to figure 2 where the block diagram corre- 
sponding to the general filter problem is given. It should be noted that 
the notation here corresponds to that in reference 4 . In the general 
situation there are usually certain elements that are given and may not 
be altered in the design of the system. These elements are denoted by 
the transfer function Ef In the figure . In particular cases the fixed 
elements might represent missile aerodynamics, servo motors, and so forth, 
depending on the application. To describe the limiting or saturating 
quantities It is necessary to express them in terms of the fixed elements 
and either the input or output of the system. Which choice is made de- 
pends on the location of the fixed elements - if the output element Is 
fixed, as is usually the case, it is desirable to express the saturating 
quantities in terms of the output. Figure 2 illustrates a situation 
wherein the output element is fixed. At this point it is necessary to 
know the saturating quantities. For this purpose we introduce a ficti- _ 
tious situation. Since in the general case the saturating quantities are 
different from the input to the fixed network, the input to Ef in 
figure 2 is imagined to be fed through certain transfer functions H sn 
which act as recording elements and whose outputs are the saturating 
quantities 0 sn . It should be pointed out that according to this formu- 
lation of the problem the outputs of the saturating elements do not feed 
Into the system but are used merely for the purpose of evaluation. From 
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the above discussion, then, it is clear that by specifying both the fixed 
and saturating networks, the saturating quantities can be related to the 
output. The remainder of the system, indicated in figure 2 by the compen- 
sating network H co is unknown and must be determined according to a 
specified criterion of merit. 


By way of comparison, the original Wiener theory considers only the 
problem of choosing the over-all transfer function so as to result in 
minimum mean- square error, 0 e , between the actual output, 6 0 > and the 
desired output, 0^* The modified theory, however, considers the minimi- 
zation of this same mean-square error with a side restriction on the 
available range of operation of the limiting quantities; for mathematical 
reasons it is convenient to consider this restriction to be in the form 
of a mean-square limitation. According to the method of Lagrange (ref. 5) 
this means that the quantity to be minimized is of the form 


©e 


Ef 



n=i 


( 3 ) 


where the p^'s represent Lagrangian multipliers which must be chosen 
properly to restrict the operation of the saturable quantities 0 sn to 
very nearly the linear range. 


In most cases it is possible to control saturation by an appropriate 
choice of the mean-square value of the saturating quantities 0 sn . In the 
simple case when a quantity has a Gaussian distribution of amplitudes the 
relation between the mean-square value of this quantity and the probability 
of its limiting is well known and particularly simple. The distribution of 
the saturating quantity In the actual case may not be Gaussian for either 
of two reasons: First the input may not be Gaussian or, second, there may 

be several saturating quantities which interact. On the other hand, even 
in these cases the distribution is frequently near enough to Gaussian that 
saturation can be readily controlled by an appropriate choice of mean- 
square value. 


The derivation of the solution for the optimum compensating network 
which minimizes the expression ( 3 ) 1b beyond the scope of this report, 
but as shown in reference ^ an expression for this network is 


Hc ° (ia ° 2*A+(m) X 


where 


-Iwt 


N 


£ 


(ia)$ ld (a)e lat 
A" (a) 


da dt 


> w 


A(co) = ^(iaOHfCiw) p n H sn (im)H SI1 (ito)ll' ii (a)) 

n=i 



8 


MCA TN 4278 


In this equation the $ quantities represent certain combinations of 
target motion and noise spectral densities (see Appendix A) , Eg refers 
to the fixed network transfer function, Es n refers to the transfer func- 
tion of the nth saturating component, and represents the correspond- 

ing Lagrangian multiplier. Of these quantities, all will normally he 
known with the exception of the . This means that only the form of 
H C o is known from this equation since the numerical constants depend on 
the values of pQ. The p n , however, are implicitly related (ref. k) to 
the mean-square restrictions on 0 sn by the following equation 


0 sn 2 = f H co (iw)H co (ia))H gn (ia))H sn (iu) 0 ii (a))da) ( 5 ) 

^ — 00 


Here it is to be noted that the integrand in equation (5) is dependent 
on the pi,p 2 , . . . ,p^ by virtue of the dependence of H c0 on these quan- 
tities (eq. (k)). It can be shown that in many cases of interest the Pn*s 
are uniquely determined by the equations (k) and (5) when a suitable mean- 

square value, 6 >sn 2 > is chosen. This quantity is to be chosen so as to 
limit the probability of saturation to a small value according to a cri- 
terion given later. In most practical applications the complexity of the 
operations involved in solving these equations does not permit the general 
solution to be obtained explicitly in terms of the input characteristics 
and fixed networks. In spite of this a great deal can be learned from 
certain numerical solutions, as will be shown. 


Application of Modified Wiener Theory to the 
Beam-Rider Guidance System 


The validity of the application of the modified Wiener theory to the 
beam-rider guidance problem depends on whether or not several restrictions 
inherent in the theory are met. In general, they are no stricter than are 
those of the unmodified Wiener theory. First, the input quantities, target 
motion and noise, must be stationary random series (see ref. 6 for a 
detailed definition) and defined by corresponding power spectra which are 
continuous. It is generally believed that displacements at the target are 
approximately stationary random series (ref. 3)* Since the beam-rider 
system operates from these displacements, the inputs to the beam-rider 
system are also stationary random series. Second, the transfer function 
of the system must have constant coefficients. In general, however, the 
kinematic loop of guidance systems involves a time -variable range factor. 

In particular, for the beam-rider system the time -variable factor Is the 
ratio of the launche r- to -mi s s ile and the launcher- to -target ranges. On 
the other hand, because the miss is determined primarily by what happens 
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near the end of flight during which the variation in this ratio is small, 
it is reasonable to assume that the requirement of constant coefficients 
is approximately met. Third, the transfer functions of the fixed and 
saturating components must he known. 

Simplifications . - There are many quantities which may saturate in the 
beam-rider system; the most important are (l) control motion, (2) rate of 
control motion, and ( 3 ) radar receiver voltages. To apply the theory to 
all these nonlinearities simultaneously would be a tremendous task. How- 
ever, there are indications that only one limiting quantity need be con- 
sidered in the present problem. 

One indication is that of the three types of limiting not all are of 
equal severity so that one type tends to predominate. Intuitively, it 
might be felt that control-motion limiting Is the most serious type and 
that placing a restriction on the control motion will satisfactorily reduce 
the other types of limiting. Of course the validity of this presumption 
would require verification. 

Another indication is that of the three types of limiting one is more 
fundamental in that it is more difficult to remove than the others . Xt 
will be noted that items (2) and ( 3 ) above are both control-system limita- 
tions while (l) is essentially an aerodynamic or mechanical limitation. 
Since control systems are relatively versatile and can be designed to cover 
a wide range of characteristics, it appears that control-motion limiting 
is the most fundamental to the problem. 

From these indications it appears reasonable to consider only control- 
motion limiting. The validity of the assumption that reducing this one 
type of limiting will reduce the other types sufficiently will then be 
verified. In terms of figure 2 this simplification means that only one 
saturating-component transfer function, H s , is involved. If the fixed 
network is chosen to be the aerodynamic transfer f miction so that the 
control motion is the input to the fixed network, then Hg equalB unity. 
Thus, a block diagram more suited to the specific case of the beam-rider 
system would appear as in figure 3* It will be noted in figure 3 that the 
input quantity called Q± in figure 2 has been split into the two parts 
which exist in the actual case, target motion y T and noise y-^. Like- 
wise the desired output of figure 2 is represented in figure 3 "by the 

target motion y^. The box Hc 0 , then, is that portion of the guidance 
system which is to be chosen to minimize the expression ( 3 ) • It may be 
noted that in the form given the diagram does not bear a direct resemblance 
to the form of an actual guidance system since the latter will involve 
feedback loops. The form shown, however, typifies that of the general fil- 
ter problem and is better suited to calculation. The solution of this 
problem can be readily converted to that of the actual beam-rider system. 
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A few additional remarks are appropriate concerning the two quanti- 
ties <t id and s>j_i in equation (4) . If 4 >i£ is defined as the autocor- 
relation function of the input in figure 2, then 

c Pii( T ) = lim bn f 0i("t)0i(t + T)dt 

T+co 


In terms of the input of figure 3> this becomes 


cpii(T) = lim 

T-»oo '-'_rp 


_1_ P T 


typ^) + yir('t)][y T (t + r) + y N (t + T)]dt 


= 9 t (t) + cp n (t) + <ftj, N (T) + ^(t) 


It Is normally assumed that there is no correlation between target motion 
and noise so that ^^(t) and are zero. It follows that the Fourier 

transfom of cp^^(t) is 

®ii(“) = ®p(^) -&■ 

Since in the present case the target motion and the desired output are 
identical it is easy to show in a similar fashion that 

$>id(“) = $p(w) 


General solution with control motion restriction .- By virtue of the 
simplifications discussed in the previous paragraphs the application of 
the theory to a beam-rider-type guidance system becomes feasible. For the 
case Illustrated by figure S> the optimum compensating network, H co , is 
given by 


Heo(i“) = 


23TA + (w) 


/ 


-ittfh r°° H f (ia)® T (a)e 

5 jL A =r- 


iat 


da dt 


where 


A (w) = Hf (ico)H f (ico) + pj$ii(w) 


( 6 ) 


In the above equation the quantities <1^ and are known from the tar- 

get motion and noise characteristics as previously discussed. The transfer 
function Hf of the missile aerodynamics was chosen for this study to 
represent that for a typical variable-incidence missile. The significant 
transfer functions for this missile are given by the following conventional 
equations (ref. 7) 
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t _ 1 %s + 1 

E Ta s(T a k s g + 2^ a T a s + 1) 

= Ta ^b 2 s g + 2^ y + i 
t Is 5 " s(%s + i) 

from which the fixed network "becomes 

= yM_ i gb a s 2 + g ^ 8 + 1 

f 5 T^ s 2 (T a 2 s 2 + 2£ a T a s + l) 


(7) 

( 8 ) 


(9) 


Table I summarizes the parameters for this particular missile at a given 
flight condition. It should "be pointed out for later comparison purposes 
that these aerodynamics correspond to those used in the previous noise 
study, reference 3- The remaining parameter p is not known so that 
equation (6) can only "be used to give the general form of H C o* To 
illustrate this, the form of Hco is derived in Appendix A where it is 
shown that 


s 2 (Ta g s g +2^ ai TgS+l) (T a g s g +2£ a T a s+l) 

(Tps+l) (T 7 2 s 2 +2£ 7 T 7 s+1) (T u 2 s 2 +2^ u T u s+ 1) (T/^-tfi^y+l) 

( 10 ) 


Many of the factors involved in this equation are functions of the unknown 
parameter p. The exceptions are (l) the aerodynamic factors, that is, 
the gain T B 2 and the second quadratic in the numerator, and (2) the terms 
due to target motion and noise characteristics, that is, the first two 
denominator factors . The remainder of the terms can only he evaluated 
after the proper value of p has "been established. This value can he 
determined from the restriction placed on the mean-square control motion. 
The analog of equation ( 5 ) in the present case is 


p 00 

B 2 = / H C o(iw)H co(iw) [Oip(u) + $jj(co)]dci> 

m 


(H) 


In theory, equations (10) and (ll) can he solved for the compensating 
network, H co , corresponding to an arbitrary mean-square control motion, 

B . In this application, however, because of the complexity of the fixed 
network and the input quantities, the computation is unwieldy. Alter- 
natively the following procedure will he used here: First calculate the 

function Hco from equation (10) for a selected value of p. After the 
result is inserted in equation (ll) , evaluation of the integral gives a 

— g 

value of S . By repetition of this process for a few judiciously 



12 


NACA TN 4278 


selected values of p it is possible to determine the value of p 
corresponding to the desired mean-square restriction g 2 ”. 

As mentioned previously, H C o does not exist, as such, in the beam- 
rider guidance system. The relationship between the filter problem being 
solved and the actual beam-rider system is illustrated by a comparison of 
figures 3 and 4. From figure 3 it is apparent that if the entire guidance- 
system transfer function is represented by Yq, then 


?o ~ S co Hf (12) 

Thus the solution for the compensating network H C o determines the desired 
over-all transfer function Yo for the actual form of beam-rider system 
shown in figure 4. 

Unfortunately, the theory does not give the performance of this opti- 
mum system directly in terms of the resultant error. The minimum error can 
be found, however, from the following relationship (ref. 7) 

e 2 = e T 2 + e-fl 2 


|l - Y 0 (iw) j 2 f^(w)da) + / jY 0 (iu>) l^jjdw (13) 

J-co 

Effects of control-deflection restriction. - The procedure discussed 
in the above paragraphs has been carried out for a range of values of the 
Lagrangian multiplier p. The results of this study are presented in 
figure 5 where the errors due to target motion and noise, the total error, 
and the control deflection are given as functions of p. These curves 
illustrate the nature of the restriction of the control deflection. Large 
values of p here correspond to a large restriction and therefore small 
control deflections. It is apparent from the expression ( 3 ) that as p 
decreases the restriction on control deflection is reduced until at p = 0 
the Wiener case is reached. By cross -plotting these curves so as to elimi- 
nate the variable p which has no significance in itself, the curve of 
figure 6 is obtained. This curve is fundamental to the problem in that it 
illustrates the dependence of error on available surface deflection. Here 
the minimum obtainable error according to Wiener theory is indicated by 
the horizontal line off to the right. The curve is shown dotted in this 
region in order to Indicate that these values of error would require unob- 
tainable-control deflections. 

The importance of the curve of figure 6 is that it represents the theo- 
retical lower limit of error corresponding to any restriction on root-mean- 
square control motion. Each point on the curve would be achieved by a dif- 
ferent guidance -system transfer function. The interesting feature of this 
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curve is that as the surface deflection is reduced from large values, the 
minimum error increases relatively slowly over an appreciable range of 
obtainable control deflections. To determine the best operating point on 
the curve it is necessary to specify the degree to which limiting -will be 
allowed to occur. Since in this linear case the control deflection is 
Gaussian, an rms control deflection of half the maximum allowable deflec- 
tion would then restrict it to within these limits about 95 percent of 
the time (ref. 8). For the specific missile being used in this study the 
maximum control deflection is about 15 which would put the desired opera- 
ting point at 7*5 as indicated in figure 6. It can be seen from figure 6 
that for this operating point the effect of this restriction is -an increase 
in error of about six feet above that of the Wiener theory. It is a some- 
what surprising result that the increase in error is so small for such a 
great reduction in surface deflection. The significance of this is, then, 
that the optimum Wiener result can be approached closely even with the 
operation confined to a realistic linear range . 

Characteristics of the optimum transfer functions .- It is of interest 
to consider in greater detail the transfer functions corresponding to the 
desired operating point in figure 6. This operating point corresponds 
very closely to a value of p = 10 4 (see fig. 5 ) which will therefore be 
used as an example. In Appendix A the transfer functions of the compen- 
sating network, closed- loop system, and open- loop system have been derived. 
The results are given below; the numerical values of the parameters are 
summarized in table II. 


T a 2 S 2 +2£ a T a s+l 

1 2 
T 2 

s 2 (T a 2 s £ +2£ a T a s+l) 

(Tps+l) (T 7 2 s 2 +2^ 7 T 7 s+1) 

T u 2 s 2 +2C u T u s+ 1 

T v 2 s 2 42£ v T v s+1 

(14) 

T a 2 s 2 +2£ a T a s+ 1 

1 

T b 2 s 2 +2^T b s+l 


(Tps+l) (T 7 2 s 2 +2£ 7 T 7 s+1) 

T u 2 s 2 +2C u T tl s+l 

T v 2 s 2 +2£ v T v s+l 

(15) 

T a 2 s 2 +2£ a T a s+l 

1 

T^ 2 s 2 + 2^T b s+l 


^ S 2 (T A S+1) 

*- JL 

T x 2 s 2 +2 £ x T x s +1 

Ty 2 S 2 +2 ^yTyS +1 

(16) 

j 


V ' 



Due to target 
maneuver and noise 
(Wiener theory) 

Due to control 
restriction 

Due to 

aero dynami c s 
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It is illuminating to trace to their origin the various factors com- 
prising these optimum transfer functions . Analysis shows that terms due 
to the basic Wiener theory also appear in the modified theory. In addi- 
tion, other terms due to the aerodynamics and the control-deflection 
restriction also appear in the solution. This situation is illustrated 
by the grouping of terms in equations (l 4 ) , (15) , and (l6) . It can be 
said that the terms which arise from basic Wiener theory correspond in 
form to the first group of terms. However, certain constants are some- 
what altered, as can be seen from table II by a comparison of the param- 
eters given for the two theories for identical target motion and noise 
characteristics. The second group of terms in the equations indicates 
that the optimum form of the filtering term due to control-motion restric- 
tion is quadratic. 2 The last group of terms in the transfer functions can 
be traced to the aerodynamics. In fact, the numerators are identical with 
terms in the aerodynamic transfer function. 

The transfer functions are illustrated in figures 7 (a ) > 7 (h), and 
7(c); these curves represent the optimum transfer functions in the pres- 
ence of the control-motion restriction. For comparison purposes the cor- 
responding transfer functions obtained from the Wiener theory are also 
plotted on these figures. It can be seen that, in general, the transfer 
functions for these two cases are similar at the very low frequencies and 
separate to a greater extent as the frequency is increased. The difference 
is primarily due to filtering terms which are required to satisfy the 
restriction on control motion. 

The added complexity has certain effects on the control motion and 
error. These effects can best be illustrated by means of the correspond- 
ing power spectra, that is, by the integrands in equations (ll) and (13) . 
These spectra have been plotted in figures 8 and 9, and for comparison the 
corresponding curves for the Wiener theory are shown. It can be seen from 
figure 8 that the spectrum of the control deflection for the more complex 
system attenuates far more rapidly than does that for the Wiener theory. 

It is as a result of this fact that the modified system does not limit. 

As for the error, it is seen from figure 9 that over the important fre- 
quency range the error due to noise is not altered too seriously so that 
most of the increase in the total error is due to an increase in the tar- 
get motion component. 

Effect of noise magnitude on performance. - The preceding results were 
determined on the basis of a design noise magnitude which lay between the 
expected extremes. The noise magnitude may vary because of targets of 
different size, aspect, and turning rate. For this reason the error per- 
formance for the optimum transfer function has been determined and Is 

3 It is interesting to note that from calculations not presented here 
it has been observed that the damping ratios of these quadratics always 
lie between 0.7 and 0.8 over the entire range_of control motion presented 
in figure 6. 
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shown In figure 10 (curve A) as a function of the zero frequency noise 
magnitude, IT. The error performance shown is opt imum only at the design 
value of noise. Shown for comparison is curve B, ■which represents the 
error performance derived from the Wiener theory, corresponding to the 
same design condition (see ref. 3 ) « The difference "between these two 
curves, then, is the increase in the error as a result of the ms restric- 
tion on the control deflection. 

As a further comparison two additional curves obtained from refer- 
ence 3 are shown. First, curve C illustrates the performance of a system 
(with certain saturable elements) which was optimized for the fastest pos- 
sible response to a step input in the absence of noise. The performance 
can be seen to be significantly poorer than the optimum performance given 
as curve A. Second, the best result that could be obtained in reference 3 
is shown as curve D in the figure. The system corresponding to this curve 
was not optimum because the linear Wiener theory used in this approach was 
not capable of evaluating all of the various filtering terms discussed in 
preceding paragraphs. For example, terms due to Wiener theory were not 
altered, and terms due to aerodynamics were not used. Terms due to con- 
trol restriction were included but the performance was limited by the 
arbitrarily assumed form. Thus the system was not optimum and limiting 
occurred; the performance deteriorated accordingly. 

System synthesis .- To achieve the results which have been presented 
it Is necessary to design the guidance system shown In figure 4 to have 
the optimum transfer function Y 0 . This can be done in many different 
ways since the actual beam-rider system consists of two distinct parts: 
the tracking radar and the missile-control system. For example, the opti- 
mum transfer function Y 0 might be split up and apportioned between these 
two parts of the guidance system. Or, on the other hand, the missile- 
control system alone might be designed to approximate the optimum Yo in 
which case the tracking radar should be designed to have a relatively fast 
response. In a similar manner the tracking radar could be optimized, in 
which case the fast response should be designed into the missile-control 
system. 

To investigate all these possibilities Is beyond the Intended scope 
of this report. As an example, however, of one possible design the choice 
of missile-control-system filtering was investigated. For this case it is 
desired to design the missile -control system to match the optimum transfer 
function Yo* This procedure is not unique; the analysis used herein Is 
based primarily on cut-and-try procedures. It has been found that the 
desired system could be synthesized in a conventional fashion illustrated 
in figure 11. To show how this system is capable of supplying the optimum 
filtering the open- loop transfer function has been derived in Appendix B 
where It Is shown that with certain assumptions 
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= ZM = ggfcgV (T 10 s + 1 )(T 4 s + l)_(^b 2 s g + 2^8 + 1) 

4 £ T d + k s k s s 2 (ThS + 1)(T 5 b + l) (as 3 + bs 2 + es + l) 

(17) 

Comparison with the optimum p D in equation (l 6 ) shows that certain dif- 
ferences in form exist. However, equation (l 6 ) can he expressed in another 
form by factoring approximately certain quadratic terms which have large 
damping ratios. Equation (l 6 ) can he written then as 

(T a s -i- l)(TqB + l)(Tb g s g + gg^B + 1) 

° S g (T^S + l) (T X S + l)[(T x S + l) (Ty 2 S 2 + 2£yTyS + l)] 

(18) 

It is now possible to choose the parameters in equation (17) to match those 
of equation (l 8 ); the specific values are given in table III. 

The remainder of the system consists of the tracking radar whose 
response should he made relatively fast for this design approach. Further, 
from the plots presented in figures 8 and 9 it is apparent that the track- 
ing radar should he designed so as not to alter appreciably the spectral 
distribution of error or control motions . Thus a break point or natural 
frequency of several cycles per second would be satisfactory. 

Comparisons and other considerations.- In a KEA.C simulation of this 
system its operation was examined in further detail. Of particular impor- 
tance in this examination was the control deflection and the associated 
servo energy. It was found that the control deflections were effectively 
restrained- to within the linear range and as a result the servo energy was 
greatly reduced. For a given servo the power expended is proportional to 
the time average of the sum of the absolute displacements of control 
motion. Thus the average servo power over a sufficient time interval of 
duration t can be determined from 


average servo power ~ — — — - 

"0 

For this system there is more than a two -thirds reduction in required servo 
energy over that for the system with smallest servo energy requirement dis- 
cussed in reference 3- The saving in servo energy becomes even greater 
when compared to systems for which the error performance becomes progres- 
sively worse than the optimum. The control-deflection time history 
obtained from the REAC was used to show that the required control rates 
are easily attainable since they rarely exceeded 3 to 3.5 radians/sec. 

This system has other virtues. For example, there is also a reduction 
of voltages within the circuit to reasonable and easily obtainable values. 
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From this simulation work, then, two facts are apparent. First, it 
is clear that the validity of applying only a control-motion restriction 
in the analysis has "been verified, since the control rate and system vol- 
tages did not saturate. Second, it is apparent that in comparison with 
systems of reference 3 in. which limiting occurred, the present system not 
only produced a smaller error "but was accompanied by other desirable 
effects such as decreased servo energy and circuit voltage requirements. 

The missile-control system discussed in the above paragraphs is 
intended only as an example designed to achieve optimum performance 
against noise. It is clear that since the system operation is essen- 
tially linear, all designs which approximate the desired Yo would give 
approximately the same performance. However, it is often necessary to 
satisfy other requirements which are related to the specific design objec- 
tives of the system and which may favor certain designs or system con- 
figurations. Because of the freedom in choice in distributing the opti- 
mum filtering in the system, the optimum characteristics specified by the 
theory do not, in general, place an inherent limitation on satisfying other 
requirements. For example, requirements on transient response during the 
beam entry phase may dictate a design in which more of the optimum filter- 
ing is located in the tracking radar. Similarly, requirements of simpli- 
city on certain parts of the system may dictate yet a different design. 
Because these and similar requirements are related to the detailed design 
objectives, an investigation of such factors is beyond the scope of this 
report . 


CONCLUDING REMARKS 


This study has considered the application of Newton’s modification 
of the Wiener filter theory to the choice of optimum transfer function 
for the beam-rider guidance system. The analysis has been applied to a 
typical variable-incidence missile at a given flight condition attacking 
a maneuvering target in the presence of glint noise. By minimizing the 
mean-square error with a side restriction on the mean-square value of 
certain quantities, limiting effects, of primary Importance in missile 
guidance, are largely eliminated. 

Although many limiting-type nonlinearities are of importance, this 
study has shown that the critical nonlinearity is due to control-deflection 
limiting, since a restriction applied to the deflection sufficient to con- 
fine control motions to within the linear range also satisfactorily pre- 
vented other types of limiting. Results have been given which illustrate 
the effect of the control-motion restriction on the minimum rms error. 

Zero restriction corresponds to the Wiener theory whose results are 
unattainable since impossibly large control motions are demanded. How- 
ever, as the surface deflection is reduced by means of an increasing 
restriction, the minimum error increases slowly until at a realistic 
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value of the control deflection the error has increased only a few feet. 
Thus it is clear that the optimum performance predicted "by the Wiener 
theory can still he approached closely hy confining operation of the sys- 
tem to a realistic linear range. 

The transfer functions of the optimum guidance system have been given 
and have been compared to the corresponding transfer functions from the 
Wiener theory. Differences were shown to be the result of added filtering 
terms demanded by the restriction on control motion. The various filtering 
factors can be traced to their origin. It was shown that terms due to 
basic Wiener theory appear in the solution, and that the additional terms 
can be associated with the control-motion restriction and the missile aero- 
dynamics. 

For the optimum system, performance against variable noise magnitude 
was given and was compared to Wiener theory. As an example of one possible 
design which achieves these results, a missile-control system was synthe- 
sized to match the optimum transfer function. In a REAC simulation the 
operation of this sytem was examined in further detail. It was found that 
the magnitude of circuit voltages, and control -motion rates are greatly 
reduced in comparison with systems in which limiting occurs . Similarly, 
the required servo energies are appreciably reduced, resulting in about 
a two -thirds saving. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif . , 


for Aeronautics 
May 11, 1955 
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APPENDIX A 


DETERMINATION OF TEE OPTIMUM MODIFIED WIENER FILTER 


The solution for the optimum compensating network H co illustrated 
in figure 3 is given hy equation (6) in the text as follows: 


H co (iw) 

where 

A(oj) 

In this equation 

®ii = ®T + (A2) 

where $rp and Ojj- represent the spectral density of the target motion and 
noise, respectively. The quantities A + and A” are defined as the factors 
of A with poles and zeros in the upper and lower half -planes, respec- 
tively. Thus 

A = A + A“ (A3) 

For the case to he considered here the target motion and noise will he 
taken as 


H-f (ia) $T( a ) e 


. - .. r e -i“tT‘ 

2jtA + (io) J 0 A (a) 

Hf (iw)Hf(icj) + pj$ii(aj) 


iat 


da dt 


> (Al) 


is defined as 


kam 

= jto^u 2 + k 2 ) ^ 

% = N (A5) 

It will he noted that the noise spectrum has heen approximated here hy a 

constant in order to reduce the complexity of the calculations. As shown 

in reference 3 , this approximation is valid because the noise defined hy 
equation (l) Is essentially flat compared to the passhand of the optimum 
system. It should also he pointed out that use of equation (Alp) leads to 
certain mathematical difficulties which occur in the process of evaluating 
the right-hand side of equation (Al) , because the theory requires that the 
poles of <&p not he located on the real axis . To avoid these difficulties 
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it is necessary to modify the target-motion spectral density to the 
following 


$rp = 


karp c 


irCw 2 + ri^) (cj 2 + t \ z z )( u £ + k 2 ) 


(AS) 


where T) x and r\ £ are any small real numbers. The solution to equation 
(Al) will then "be a function of iw, and t) z . The desired answer is 
obtained by taking lim Y 0 (ico,^,^) as and r> -> 0. However, it can 

be shown (see ref. 3 ) that the same answer can be obtained more simply by 
taking t) and r\ equal to zero as would be obtained from the rigorous 
process described above . 

Now' by combining equations (A2), (Ak) , and (A5) , 


$ii(a) = ^(a) + %(a) 

a 6 + k 2 a 4 + (karji 2 /irN) 
a 4 (a 2 + k 2 ) 


(a - a P i)(a - ap g )(a - aps) (a + cc px ) (a + apg) (a + aps) 
a 4 (a + ik)(a - ik) 


[p(cx)3[-p(-a)3 

a 4 (a + ik) (a - ik) 


(A7) 


where the apm represent roots in the upper half -plane. The polynomial 
p(a) is 

3 

p(a) = n (a - cc^) 

m=i 

= a 3 + b e a 2 + bia + bo (A8) 

where 

b 2 = “(api + ap 2 + aps) 
bi = apxOpg + ap X ap3 + ap 2 aps 
bo = "Opiap 2 c£p3 


> 


(A9) 
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In a similar manner 


-p(-a) = a 3 - b 2 a 2 + bid. - "bo 


(A10) 


An expression to be used later can be obtained from these preceding 
equations as follows 
ka 2 

a 6 + k 2 a 4 - 1 — = (a 3 + b 2 a 2 + bia + bo) (a. 3 - b 2 a 2 + b x a - bo) 


from which 


2 _ 


“bo = 


jrN 


(All) 


The quantity Hf> in equation (Al) represents the fixed network or 
aerodynamics and its general form (ref. 7) is given by 


„ f. S _ *M _ 1 [(1 - T b a a g ) + ite^yQ] 

j & Ts 2 " a 2 [ (l - T a 2 a 2 ) + i(2£ a T a a)] 


(A12) 


and 


tt ' 7 7— N _ 1 [(1 - T^a 2 ) - i(2e- h T- h cx)3 

Hf ( 1 ) Tp - a 2 [ (l - T a 2 ct 2 ) - i(2? a T a a)] ^ A13 ' 


Thus 


5-Tr-^w t* \ - 1 [Tb 4 a 4 + (^gh^b 2 - gTb 2 )^ 2 + 13 
f( ) f( ) T s 4 a 4 [T a 4 a 4 + (l^ 2 ^ 2 - 2T £ 2 )a 2 + 1] 


(AUO 


Accordingly it can be shown that in equation (Al) 


Hf (ia)Hf (ia) + p 


_ pT s 4 Ta 4 a a +pT s 4 (^-g a 2 T a 2 -2T a 2 )a, 6 +(T b it +pT s 4 )g‘ t +(^g-b 2 Tb 2 -2Tb 2 )a 2 +l 
T s 4 oj 4 [ T a 4 a 4 +(l|.£ a 2 T a 2 -2T a 2 ) a 2 +l] 

II (a-oqn) II (aKCq n ) 

= nT 4 n=i n=i 

P a a 4 [ T n 4 a 4 +( 4 f 2 T 2 -2T 2 ) a 2 +l] 


= pT a 4 


q(a)q(-a) 


a 4 [ T a 4 a 4 + ( 4 £ a 2 T a 2 -2T a 2 ) a 2 +l ] 


(A15) 
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where the Oq n represent roots in the upper half -plane. The polynomial 
q(a) above is 

4 - 

q(a) = II (a - aqn) 
n=i ^ 

= a 4 + c 3 a s + c 2 a 2 + o x a + c Q (Al 6 ) 

where 

c s = - (ccqi + aq 2 + aqs + aq 4 ) 

c 2 = aq 1 aq 2 + aq 1 aq S + aq 1 aq 4 + aq 2 aq 3 + <xq 2 aq 4 + aq 3 aq 4 

> (A17) 

c i = " (<XqiCCq 2 Oq3 + ocqiOq 2 aq 4 + Oq 1 aq 3 0 Lq 4 + ctq 2 ctq 3 aq 4 ) 
c Q = otqiaq 2 aq 3 aq 4 
In a similar manner, 

q(-a) = a 4 - 030 ,® 4 - c 2 a 2 - Ca.a + c 0 (A18) 

An expression to he used later can he obtained from these preceding equa- 
tions as follows 

pT s 4 T a 4 a s +pT s 4 (4? a 2 T a 2 -2T a 2 )a 6 +(T b 4 +pTg 4 )a 4 +(4^ 2 T- b 2 -2T b 2 )a 2 +l 

= pT g 4 T a 4 q(a)q(-a) 

= pTs 4 Ta 4 (a? 1 +c 3 a a +c 2 a s +c 1 a+co) (a* -c 3 a 3 +c 2 a 2 -c 1 a+Co) 


from which 


Co 2 = 


pW 


(A19) 


The second equation in (Al) can now be found by combining equations (AT) 
and (A15) 


A(a) at pT a 4 N 


q(a)q(-a)p(a) [-p(-a) ] 


a. s [T a 4 a 4 + (4^ a 2 T 2 - + 1 ] (a + ik) (a - ik) 


(A20) 
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This expression can he split into two factors with poles and zeros in the 
upper and lower half -planes . From the preceding definitions q(a) and 
p(a) both have their zeros in the upper half -plane. Thus 


A + (q) = pTa 4 N 


q(a)p(a) 


a 4 (l - T a s a 2 + i2^ a T a a) (a - ik) 


(A21) 


A" (a) = 


q(-a)[-p(-a)] 


a 4 (l - T a 2 a 2 - i2^ a T a a) (a + ik) 


(A22) 


Now, let us consider the evaluation of the first integral in equa- 
tion (Al) . This integral can he evaluated by combining equations (A4) , 
(A13) , and (A22) as follows : 



Ef (a) $p(a)e j ' a ‘^ 
A”(q) 


dq 


*°T e r” a - - ia^a)e iat 

Ts £ tt J_ m a? (a - ik)q(-a) [ -p(-a) ] 

= _ r“ (1 - yq 2 - i2^q)e lat ^ 

T s k -00 a S( a _ i^) n (q + qqn) II (q + qpm) 

n=i m=i 

kam 2 r°° 

= - f(q)d q (A23) 

T g It <J_oo 

where f(q) is defined as the integrand. By considering q to he a com- 
plex variable, the integral can he evaluated hy a contour integration indi- 
cated hy the following sketch: 
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There are two poles within this contour, a second-order pole at the origin 
and a simple pole at a = ik, as indicated. The residues are 

Res(o) = lim [a 2 f(a)] 

' ' a*o da v J 

t khpCpRi^Tfo + ik(bpc x + b^ep) + bpcp 
khpCp k 2 h 0 2 c 0 s 


= + 7 X 

kb 0 c 0 


(A24) 


Res(ik) = lim (a - ik)f(a) 
a-*ik 

= (T^k 2 + gg>,Tbk + l)e~ kt 
-k 2 q( -ik) [ -p( -ik) ] 


s r 2 e" kt (A25) 

Now, hy means of the residue theory, the value of the integral is given 

t>y 


2ikam a 

Ii ^ — [Res(o) + Res(ik)] 

Tc 


Thus from equations (A24) and (A25) , 


I, = 

T S 2 


^ + i^ +72e " kt 


0 


(A26) 


The second integral in equation (Al) is denoted hy I 2 and is merely 
the Fourier transform of I 1# Hence 


-f 


I 1 e" i “ t dt 


2ikap £ 

T s 2 


y p -icot . — £ — e -±wt , y e -(k+: 

r i e + kh 0 c 0 e + 7 s e 


iGJ ) t Jat 


(A27) 
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Thus it can "be shown that 

2 ika T £ ( 73 . + 7 g)b 0 c 0 (ito) 2 + [ky 1 h 0 c 0 + (l/k)3ia) 4 - 1 
2 “ " T E 2 b 0 c 0 (iu>) 2 (k + iu) 


T CL 2 (iw') 2 + 2£ a T a (lu) + 1 
T s 2 b 0 c 0 (iw) 2 (k + iw) 


(A28) 


It should he noted that since the T] , s were assumed to he zero, questions 
concerning the existence of the above integral arise. However, as indi- 
cated in reference 3* "when the analysis is made without this assumption, 
the integral in equation (A27) does exist, and the limit of this integral 
as the ti’s approach zero becomes precisely equation (A28) . 

The coefficients can he expressed in a more convenient form by elimi- 
nating the intermediate parameters y x and y 2 by means of definitions 
given in equations (A24) and (A25) . The following results are then 
obtained: 

2 ?a T ct = hy^QCo + 


- - 1 ' b ° C \oCo lg ° (A29) 

Also 

T a 2 = (7i + y 2 )h 0 c 0 

m _ ‘ b o c o(^ ak2 + 2 £hTfr)K + l) _ kb 0 c 0 2^Tb + ikCbpCx + bxCp) _ 1 
k 2 q(-ik) [ -p(-ik) ] k^oCo k 2 

(A30) 

For the missile parameters of table I and the value of k given on page 28 

T b 2 k 2 + 2£ b T- b k + 1 « 1 


so that very nearly 


T a 2 


~bp c o 

k^qC-ikH-pC-ik)] 


_ 1_ 
k k 2 


+ 


(A31) 
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Now, from equation (Al) the optimum compensating network can "be 
determined as follows: 


EcoC^) 


(A32) 


2jcA + (u) 

Substituting equations (A2l) and (A28) into (A32) gives 
Hco(i w ) 

karp 2 ( m) 2 [T a 2 ( ico) 2 +2^T n (iu/) +1] [T r 2 (1 m) 2 +2 g a T a (iu) +1] 

3tNpTg 2 Ta 4 t>o c o p(w)q(a)) 

(A33) 

To reduce this expression to a more standard form it can "be seen from 
equation (Al6) that 

q(“) = n 3i (« - a qn) 

- °° nl ('4r + ^ (A3*) 

Similarly from equation (A8) 

3 

p (“) = m Si 

= h 0 S (- gS - + l') (A35) 

u m=i \ a pm / ' 

Also hy using the definitions given in equations (All) and (A19) the com- 
pensating network can he shown to reduce to 

HcoUw) = T s 2 ( la) ) 2 [ T a 2 ( i ^) 2 + gSafti(m) + l3[T a 2 (i(d) 2 + 2£ a T a (ito) + 1] 


jSx (" a pn + 0 n3i (" aqn + ^ 


(A36) 


By the substitution of thd conventional complex frequency s for iw 


Hco(s) „ Ts a ± + 1) ( V B g + gta T a° - D 


3 

n 

m=i 


Si + OnSi fe + ) 


(A37) 
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It is shown later in the numerical examples that certain terms combine to 
give quadratic factors as shown below 


Hco ( s ) - Ts‘ = 


s a (l a g s 2 4-2g fT T a s+l) (T a 2 s 2 +2^T a s+l) 


(Tps+l) (T 7 2 s 2 +2^T 7 s+l) (T u 2 s 2 +2^T u s+1) (T v 2 s 2 +2£ v T v s+1) 

(A 38) 


where 


T 0 = 

X 

apx 

T 2 
X U 

II 

EM 

1 

a P2 a pS 

2£u T u 

2£ 7 T 7 - 

^ ap2 + a,ps 
ap 2 aps 

Tv 2 



2£v^v 


1 

«Xqictq2 

gqx + gq2 
“qittqa 

1 

<Xq3aq 4 

Qqs + Qq4 
a q3 a q4 


The compensating network given in equation (A38) can he used to deter- 
mine the over-all transfer function Y 0 . From figure 2 

Y 0 (s) = H co (s)H f (s) (A39) 

Substitution of equations (A12) and (A38) into (A39) gives 


Yo(s) 


(T^ 2 s 2 +2g h T^s+l) 

(Tgs+l) (T 7 2 s 2 +2^ 7 T 7 s+l) (T u 2 s 2 +2^T u s+1) (T v 2 s 2 +2£ v T t s+1) 

(AllO) 


If Yo is considered to represent a unity feedback system, the equivalent 
open-loop transfer function is shown in the later specific example to be 
of the following form: 


M-o(s) 


Yo(s) 

1 “ Y 0 (s) 


v (T a 2 s 2 4- 2g„T„a 4- 1) (T^s 2 4- 2^8 + 1) 

^ s 2 (T^s 4 - 1)(T x 2 s 2 + 2£ x T x s+l) (Ty 2 s 2 + 2£ y TyS + l) 
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Because and £ x are close enough to unity (table II) it is possible 
to factor the corresponding terms approximately to give 

* , (V + l)(l tt Sfl)(^^2 Wt l) 

U s 2 (T x s + l)(T x S + l)[(T X S + l)(Ty 2 S 2 + 2CyTyS + 1)] 

(Ate) 

The numerical evaluation of the optimum, system has been carried out 
for the following values of target motion and noise: 

k = 0.4 ®T = 1 S 

T = 5 sec N = 15 ft 2 /radian/sec 

From this example p will be taken to be 10 3 4 * since, as shown in the text, 
this value results in a near-optimum restriction on control deflection. 

Evaluation of the numerator in equation (A7) gives 

ka ™ 2 

a 6 + k 2 a 4 +— rc = a 6 + 0.l6 a 4 + 8.8 

3 3 

= ( a “ a pm) ( a + a pm) 

where 

a px = i 1.456 

ap 2 = 1.427 exp(i O. 543 ) 

a ps = -1.427 exp(-i 0.543) 

which are a] 3 located in the upper half of the complex a plane. Then 
from equation (A9) 

b 2 - -i 2.94 
b x = -4.201 
b 0 = i 2.97 



MCA IN 4278 


29 


Similarly from equation (A15) for p = 10 4 , 

pT s 4 T a 4 a 8 +pT s 4 (4^ a 2 T a 2 -2I a 2 )a 6 +(l b 4 +pT s 4 )a 4 +(l;^ 2 I b 2 -2I b 2 )a 2 +l 

= (22.58xl0" s )a 8 -C7. I )-75xl0" 5 )a e +(6.267X10“ s )a 4 -(0.00609it)a 2 +l 

= pT s 4 T a 4 ^ (a-ctqn) (a+aqn) 


Solving for the roots gives 

ttqi = 3.593 exp(i O. 80 I) 

Oq 2 = “3*593 exp(-i 0.801) 
aqs = 12*92 exp(i 0 . 058 ) 
aq 4 = - 12.92 exp(-i 0 . 058 ) 

•which again are all in the upper-half a plane. From the definitions in 
equation (A17) , 

C 3 = -i 6 . 6 77 

c 2 = -I 87 . 6 I 
Ci = i 883.10 
Co = 2155 


Now from the values of the roots ctp^ and Oq n as well as the h an d c 
coefficients, the optimum compensating network can he calculated hy means 
of equation (A 38 ) . For this purpose it is necessary to evaluate the left- 
hand sides of equations (A29 ) , (A 31 ) , and the parameters following equa- 
tion (A38) . ■ Results of these calculations are given in table II in terms 
of damping ratios and undamped natural periods. Substitution of these 
parameters in equations (A 38 ) and (A4o) gives 



CO 

o 


Hco(s) = 7.91X10" 4 


g g (l,28 S g +1.82 B+l) (Q. 00602 B g +0 ,00830 B+l) 

(0.687 s+l) ( 0.490 s 2 40.727 s+1) ( 0.0775 s 2 + 0.401 s+l) (0.006 s 2 + 0.009 s+l) 


(A43) 


T g (s) 


(1.28 b £ +X, 82 s+l) (0.00305 s 2 +0. 00243 s+l) 


(0.687 s+l) (0.490 b 2 +0.727 s+l) (0.0775 s 2 -K).401 s+l)(0.006 s 2 +0.009 s+l) 


(a44) 


AIbo Included. In this table are the parameters of the- equivalent open-loop transfer function 
(eq. (A4l)) which Is 

(1.28 s g + 1.82 s + 1)(0.0Q305 s g + 0.00243 s+l) 

Ho(s) ■ 3-2 s a( 2<5 g + !) (0.0338 s 2 + 0.266 s + 1) (O.OO619 b 2 + 0.00807 s+l) 

(A45) 
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APPENDIX B 

MISSILE-CONTROL -SYSTEM APPROXIMATION TO 
THE OPTIMUM TRANSFER FUNCTION 


The optimum linear system is illustrated in figure 11. In the follow- 
ing section the system equations are derived. 

Using figure 11 it can he shown that 

Jfr_ _ ^s(^M s + 3-) (^3 S + l) 

% T d s(T lS + l)(T a 2 s 2 + 2£ a T a s + l)(T s s + l) + k s k s s(T M s + 1 ) (T 2 s + l) 

(Bl) 

To simplify equation (Bl) it is convenient to make T 3 = T^. This assump- 
tion is not essential hut its use leads to simpler equations . In certain 
cases where gust disturbances are serious it may he more desirable to 
choose T 3 small. With the former choice 

_ k s (T M s+ 1 ) 

% s[(T d T a 2 T 1 )s 3 +(T d T a a + T d 2^ a T a T 1 )s 2 + (T d T 1 + T d 2^ a T a +k s k s T 2 )s + (T d +k 6 k s )3 

- 

T d +k s k 3 


(%s+l) 

s (as s +bs 2 +cs+l) 


(B2) 


where 


a = _T fl T g 2 T i _ 

+ ^6^3 

t - T *( T a 2 + g Wi) 

T d + k s k 3 

c _ 2 ^a T a T d + T i T d + k s k s T g 
T d + k s k 3 


(B3) 
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Thus the entire open-loop transfer function can be written as 

_ yM _ k g k 8 V (T 10 s + 1)(T 4 S + l)(T b g s g + 2£ b T b B + l) 

^ 1 T d + k s k s b 2 (t i;l s + l)(T 5 s + l) (as 3 + bs 2 + cs + l) 

(B4) 

By comparison with equation (A42) it can be seen to be of the same form as 
the optimum transfer function. Thus the following correspondences between 
quantities in equations (b 4) and (A42) are appropriate: 

Tn <-* T^ 

T 5 «-> T x 

> (B5) 

T 10 «— » Toe 
T 4 <— > T<x 

* 

Also it is apparent that 

a T x Ty 2 

b T y 2 + T X 2^yTy > (B6) 

c «->T X + 2£yTy - I 


For these values of a, b, and c, equations (B3) can be solved for the 
remaining parameters to give v 


m = aT 'a 

1 - a2£ a 


k s k 3 = 


TdTa 2 Tx - aT d 


> 


(B7) 


Ts = cT d + ck s k 3 ~ 2 ^a T a T d - T i T d 
k 8 k 3 


J 



NACA TN 4-278 


33 


The value of k 2 k B is determined from the desired gain: 

+ kpko , 

k 2 ks <-> kp. a- - S - Q (B8) 

Since only k 2 k B and k s ks are specified, one of the three gains k 2 , k 3 , 
or ks may he chosen arbitrarily, subject to the condition that voltage 
limiting does not occur. From the above equations the system parameters 
have been determined and are tabulated in table III. 
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TABLE I.- SUMMARY OF AERODYNAMIC PARAMETERS FOR EXAMPLE MISSILE 


Parameter 

Value 

T a 

0.0775 


-0552 

Td 

2.087 

T s 2 

.0007911 

% 

.846 

Ca 

.0536 

£b 

.0220 


TABLE II.- SUMMARY OF PARAMETERS OF OPTIMUM TRANSFER FUNCTIONS 


Parameter 

Modified Wiener 
theory value 

Wiener theory 
value 

T a 

1.13 

0.925 


.805 

.765 


.687 

.687 

t 7 

• 700 

.700 


.519 

■519 


2-5 

2-5 

ky. 

3.2 

7.42 

Tu 

.278 

— 


.720 

— 


.0774 

— 


.0581 

— 

T x 

.184 

— 


.722 

— 

% 

.0787 

— 


.0513 
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TABLE III.- SUMMARY OF PARAMETERS OF OPTIMUM SYSTEM 


Parameter 

Value 

k 2 

0.0538 

t 4 

1.13 

T10 

1.13 

T11 

2-5 

t 5 

.181+ 

k S 

.0463 

Ti 

.186 

k 3 

-.674 

T 2 

• 373 

Ts 

.846 





(b) Block diagram. 

Figure 1.- Beam-rider guidance system. 
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3 .- Block diagram of filter problem for the hearn-rider system. 


Missile-control system 



Figure 4.- The heam-rider system. 






RMS error components, feet 




Lagrangian mu (ti pt ier p Lagrangian multiplier p 

(a) RMS error. (b) RMS surface deflection. 

^Figure 5 - ~ Variation of error and control motion with Lagrangian multiplier p. 
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RMS error feet 




Figure S.- Effect of rms control motion on minimum rms error. 
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Figure 7 .- Optimum transfer functions. 
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Frequency, cps 

(b) Closed -loop transfer function Y 0 . 


Figure 7*- Continued. 


Phase angle , degrees 
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(c) Open -loop transfer function JJL 0 - 


Figure 7*“ Concluded. 


Phase angle, degrees 
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Figure 8.- Spectral density of control motion 
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Figure 9.- Spectral density of error components. 
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Figure 
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Comparison of ms error 
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performance of modified Wiener theory 
; optimized for design value. 























